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ABSTRACT 

Classical Novae (CNe) are the brightest manifestation of mass transfer onto a white dwarf in a cataclysmic 
variable (CV). As such, they are probes of the mass transfer rate, M, and WD mass, Mwd, in these interact¬ 
ing binaries. Our calculations of the dependence of the CN ignition mass, Mign, on M and Mwd yields the 
recurrence times of these explosions. We show that the observed CNe orbital period distribution is consistent 
with the interrupted magnetic braking evolutionary scenario, where at orbital periods Porb > 3 hr mass transfer 
is driven by angular momentum loss via a wind from the companion star and at Porb < 3 hr by gravitational 
radiation. About 50% of CNe occur in binaries accreting at M ~ 1O“‘^M0 yr“' with Poib = 3-4 hr, with the 
remaining 50% split evenly between Porb longer (higher M) and shorter (lower M) than this. This resolution 
of the relative contribution to the CN rate from different CVs tells us that 3(9) x 10^ CVs with WD mass 
l.O(O.6)M0 are needed to produce one CN per year. In addition, one CN per year requires a CV birthrate of 
1(2) X 10“^ yr“*, and likely ejects mass into the ISM at a rate M = 3(9) x 1O“^M0 yr“'. Using the K-band 
specific CN rate measured in external galaxies, we find a CV birthrate of 2(4) x 10“'* yr“* per 1O'‘’L0 ,a:, very 
similar to the luminosity specific Type la supernova rate in elliptical galaxies. Likewise, we predict that there 
should be 60-180 CVs for every 1O*’L0 _a: in an old stellar population. The population of X-ray identified CVs 
in the globular cluster 47 Tuc is similar to this number, showing no overabundance relative to the field. 

The observed CN Po^b distribution also contains evidence for a CV population which has no period gap. 
These are likely systems with a strongly magnetic WD (Polars) in which it has been suggested that the field 
interferes with the wind of the companion, limiting angular momentum losses to those of gravitational radiation 
and eliminating the period gap. With this reduced M, Polars evolve more slowly than systems which undergo 
magnetic braking. Using a two-component steady state model of CV evolution we show that the fraction of 
CVs which are magnetic (22%) implies a birthrate of 8% relative to non-magnetic CVs, similar to the fraction 
of strongly magnetic field WDs. 

Subject headings: binaries: close—novae, cataclysmic variables- stars: statistics —white dwarfs 


1. INTRODUCTION 

A Classical Nova (CN) outburst is the result of a ther¬ 
monuclear runaway in the hydrogen-rich accreted layers on a 
white dwarf (WD ) in a mass tr ansferring cataclysmic variable 
(CV) binary (see IShar^l 19891 for a review). The number of 
CNe with mea s ured o rbital periods, Porb, has been increasing; 
iDiaz & Bruchl (Il997h made a first compariso n to the CN pe- 
riod distribution with 30 objects, and recently FWarneri (l2002li 
listed Porb for 50 CNe. These range from 1.4 hours to more 
than 16 ho urs, with most below 5 hours. Cataclysmic vari¬ 
ables (CVs: IWarneTl199^ consist of a WD primary accreting 
mass from a normal low-mass star that is filling its Roche 
lobe, and spend mo st of their life accreting at rates M k, 1 0“^- 
1O“"M0 yr“* (e.g. [Howell, Nelson, & Rappaportll2001[ her- 
after HNR). CVs are formed when the companion (of mass 
M 2 ) to a WD (made during a common envelope event) comes 
into contact with its Roche lobe as a result of angular momen¬ 
tum losses. The CV evolves toward shorter orbital periods as 
the binary loses angular momentum at a rate J and the com- 
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panion transfers mass to the WD at a rate M. The orbital 
period distribution of CN, a direct reflection of the CN fre¬ 
quency, provides a diagnostic of M in CVs over the time to 
accumulate an unstable layer, 10^“’ years. 

It is generally th ought that the paucity of CVs with 2 hr 
^ ^orb ^ 3 hr (e.g. l■ShafteT^^1992h . the so-called period gap, 
is due to switching between a high J state above 3 hours 
to a comparatively low J state below 3 hours (e.g. HNR). 
Upon such a sudden reduction in J, the requisite reduction 
in M causes the companion, which has been slightly out 
of thermal equilibrium due to the high mass loss rate, to 
contract and fall inside the Roche lobe, halting mass trans¬ 
fer. Eventually as Porb continues to decrease due to an¬ 
gular momentum loss, Roche lobe contact is reestablished. 
The simplest paradigm of CV evolution holds that J is de¬ 
termined by the WD mass Mwd, M 2 and Porb. Interrupted 
magnetic braking (IMB) has arisen as a good candidate 
for J prescriptions (see e.g. [^czynsk i&Me nkiewicz y^S^ 
iRappaiport et al.l[T90 ISpruit & Ritterlll98^Hameurv et alJ 
rr^ Ik^dUrr^ HNR). in the IMB scenario, J is ini- 
tally driven by angular momentum lost in the stellar wind 
of the secondary star (IVerbunt & ZwaarJ 119811) for Porb > 
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3 hours. The magnetic braking is interrupted (potentially 
by the loss of the radiative core), and J f alls to that car¬ 
ried away by gravitational radiation (see e.g. |Faulkn eHll97U 
iPaczyn^i & Sienkiewiczil98lURaDpaDort et al.l 19821) . IMB 

explains the period gap of the CV population and the pre¬ 
dicted jM) agrees with observat i ons of the WD surface tem¬ 
peratures dTownslev & Rildstenll2003li . This agreement is 
somewhat surprising because modern measurements of J for 
isolated low mass st ars do not appear to fo llow the law used 
in the IMB scenario dAndronov et al J2003I) . IMB also makes 
a very important but unconfirmed prediction: that CVs evolve 
across the period gap. 

For a given composition of accreted material, there are three 
independent parameters that set Mign, and thus the CN rate: 
M, 7) , and Mwn. Many previous studies (e.g. lTruran & Livi^ 
I198 ct iRitter et ^119911) focused only on the dependence of 
the CN frequency on Mwd- Here we focus on T^’s depen¬ 
dence on M. iTownslev & BildstenI d2004li (hereafter TB04) 
showed that, after extended accretion and many CN cycles, Tc 
is set by M because of accretion’s impact on the thermal state 
of the WD core during the long interval between CNe. Given 
any M(Porb) relation, TB04’s calculation of Tc enables the first 
consistent prediction of how the CN rate varies with orbital 
period. Usin g M dependent M\„„ values calculated at sev- 
eral Tc values dPrialnik & Kovetall995h . iNelson et alJ (l2004l) 
showed clearly that this knowledge of Tc is essential for a con¬ 
clusive prediction of the Porb distribution. In Section |2l we 
briefly review the work of TB04, and in Section 0 we apply 
it to calculate the CN period distribution and compare to that 
observed. We And that IMB is remarkably consistent with the 
observed CN frequency, yielding a Arm relationship between 
observed CN frequency and the CV population. Confirming 
this scenario Anally allows us to use the observed CN rate to 
infer the number of underlying CVs. This impacts searches 
for CVs in the held, as well as in globular clusters, as we note 
in Section^] 

2. CLASSICAL NOVA IGNITION MASSES 

IMB predicts that typical accretion rates near Porb = 3- 

4 hours are M ~ 1O“‘’M0 yr“', and decrease to M ~ 5 x 
1O“"M0 yr“' below 2 hours. This factor of 20 drop in 
M across the period gap requires that the M dependence of 
Mign be taken into account when explaining the orbital pe¬ 
riod distribution of CN. These two M’s lead to different ig¬ 
nition modes. For the higher M, TB04 found t hat the CN i s 
triggered by trace ^He in the accreted material (IShar ’MM, 
leading to a CNO burning runaway in a radiative layer. At 
the lower M, however, the CN outburst is triggered by the pp 
chain when the base of the accreted shell is conductive. This 
wide range of conditions means that the pr ior prescriptions of 
ignition, such as a single unique pressure llGehrz et al.lll998l 
WarneJl2002h or taking M {„„ independent of M (iRitter et al.l 
199lUDiaz&Bruchll997ll . are wholly inadequate. 

By considering the release of energy due to compression of 
the accreted layers, TB04 calculated the rate of heat loss or 
gain of the WD core. This rate varies as the accreted layer 
builds up, and by taking the average loss or gain to be zero 
TB04 found an equilibrium, Pc.eq, which the WD will ap¬ 
proach under prolonged accretion. For Porb ^ 3 hours, the 
magnetic braking phase, the WD has only marginally enough 
time, ^10^ years, to reach rc,eq- However, at these accretion 
rates, Mign is insensitive to Tc if it is lower than Pc^eq because 



Fig. 1.— Classical Nova ignition masses, Mign, for WDs which have 
reached their equilibrium core temperature, as determined by TB04. Con¬ 
tours are evenly spaced in log(Mign), labels indicate M^gn in Mq. The 
vertically hatched region indicates where steady burning of H is expected 
IPiersanti et alllTOOCI: In omo laiim. At higher M, envelope buildup and 
expansion into a giant is expected. 


the ignition occurs in a layer which is not conductively cou¬ 
pled to the core. Hence the temperature is set by the response 
of the envelope to the energy being released by compression 
(TB04). At lower M, Tc determines Mign, but the low M state 
for CVs is long-lived (> 10® years), allowing the WD core 
to reach Pn eq. Thus, for CV evolution, use of P^ eq provides 
Mign’s which are self-consistent and independent of assump¬ 
tions about the WD thermal sta te. This is in contrast to all 
prior Mi„n calculations (see e. g. lFuiimotolll982t iMacDonaldl 
ll^lPrialnik & Kovetdl^ . The Tc values found by TB04 
are <10^ K for Mwn < 1.2 Mm, low er than the lowest Tc con¬ 
sidered bv iPrialnik & Koveta ill 99.51) in their broad parameter 
survey. 

Models of donor evolution in CVs llD’ Antona & M azzitellil 

Il982tliben & Tutukovll984tlMcDermott & Taamll989l) show 

that the ^He abundance reaches mass fractions of Ashe — 
0.001 when the orbital period is 3-4 hours (the compan¬ 
ion mass is below ~ O. 4 M 0 ). This is when the star has 
lost enough mass and the convection layer has become deep 
enough that the partially-fused material in the stellar interior 
is uncovered. Precisely when this occurs and the As^e reached 
is sensitive to the donor’s initial mass, but generally does not 
exceed 0.003. We are most interested in CN rates for Port < 4 
hours, so we include a fraction Asue = 0.001 in our ignition 
models. 

The Mign results of TB04 are shown in Figure^ Compared 
to previous ignition mass calculations, TB04 found larger 
Mign at M < 1 O“'°M 0 yr“^ due to the lower Tc, and simi¬ 
lar Mign at Mign ^ 1 O“®M 0 yr“' due to offseting effects of 
the lower Tc and the inclusion of ^He leading to an earlier 
trigger. TB04 found that for the few systems which have a 
known ejected mass and Porb (giving some indication of M), 
the ejected masses are similar to Mign, implying that each CN 
ejects approximately the amount of material accreted since 
the last event. TB04’s results do not show the trivial Mign oc 
P‘*/Mwd scaling assumed for a unique ignition pressure (see 
their discussion). An ignition pressure of p = 2 x 10 '® dyne 
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Fig. 2.— The cumulative (top panel) and normal (bottom panel) distribu¬ 
tion of classical novae orbital periods. The thin solid lines show the observed 
distribution, with the 9 systems with Porb > 6 hours not shown. The steady 
state population described in the text {thick solid line), consists of an equal 
number of CVs which undergo an interrupted magnetic braking episode (IMB 
population, dot-dashed line) and which undergo no magnetic braking (NMB 
population, dotted line). This combined population reproduces the major fea¬ 
tures of the cumulative distribution for Porb ^ 4 hours. A full population syn¬ 
thesis (HNR, dashed line) including only objects which undergo a magnetic 
braking phase displays similar features, matching the observed distribution 
remarkably well. A mass fraction = 0.001 is used for all. 


cm“^ would have 10“^, 10“^ and IQT^Mq contours which are 
vertical lines at M wn = 1.3. 1.0 and0.58M©. (This ignition 
pressure is used bv iGehrz et al.l lll998h in their discussion of 
the nucleosynthetic impact of nova ejecta.) As we will see, 
most CNe have M ~ lO’^M© yr“*, and thus lower Mign than 
this prescription implies. 

Also shown in Figure Q is the region (hatched) where 
steady burning is e xpected. For 0 . 52 < Mwn/Mo) < 0.68, 
the calculations of iPiersanti et al.l il20f)0l) are shown. For 
Mwd > 0.6Mq the u pper bound is found from the core mass- 
luminosity relation of iPaczvhskil Jl970ll . and the lower bound 
is taken to be 0.4 of this d omo tnll982h . At higher M, enve¬ 
lope buildup and expansion into a giant is expected. 

3. THE CLASSICAL NOVAE ORBITAL PERIOD DISTRIBUTION 

The observed distribution of CN orbital periods dWarned 
l2002h is shown in Figure |2lby the thin solid li ne. These data 
differ only very slightly from the collection in iRitter & Kolbl 
d2003l) . We have neglected any selection effects, as they have 
not been satisfactorily quantified to date. These are complex, 
so we will not attempt a conclusive discussion, only men¬ 
tion two important ones, C N outburst brightness and remnant 
brightness. The tables in iPrialnik & Kovet j d 19951) . for the 
coldest Tc and Mwd = IM©, imply that CN outbursts for sys¬ 
tems with M ^ 10“'°M© yr“' are 3-4 times as bright as those 
for systems with M ~ 10“*^M© yr“*. However, the spatial 


distribution in the Galaxy and the effects of absorption must 
be considered in order to judge how th e observed populat ion 
reflects this brightness variation. (See IRitter et al.ll 19911 for 
some discussion.) It is likely that remnant brightness depends 
on M, and thus Porb. introducing additional selection effects in 
remnant recovery and characterization. 


3.1. Using Previous Population Synthesis Calculations 
In order to calculate the CN frequency we must first use the 
IMB prescriptions to obtain M from Porb- Since the secondary 
fills its Roche Lobe, Po^b is directly related to M 2 by its mass- 
radius relation. We enforce R 2 = 0.46a[M2/(MwD + 4f2)]*^^. 
so that Porb = 9 hr {M 2 /Mq)~^^^{R 2 /. To reproduce 
IMB for Porb > 3 hours we use the magnetic brak ing law from 
HNR (see also lRappaport, V^bunt, & Jossll98^ 


■/mb = -9.4 X 10^^ ergs 


Ml.] 

Me) 



( 1 ) 


When the binary is evolving under the influence of mag¬ 
netic braking, the M 2 is large enough that the secondary is 
out of thermal equilibrium and therefore has a larger radius 
than a main sequence star of the same mass (the effect that 
leads to the period gap). In order to account for this we 
use the M 2 -R 2 relation found by HNR using a bipolytrope 
model for the secondary and evolving under equation Q, 
R 2 /RQ = O.81(M2/M0)°-67. For Po* < 3 hours, we apply 
gravitational radiation losses given by 


7gr- 


32GQ^uj^ 

5c^ 


= -2.7 X 10^'^ ergs 



MwdM2 

M,M© 


2 



,( 2 ) 


where Q is the moment of the binary about the orbital center, 
02 = 27r/Porb, a is the orbital separation and M, = Mwd + M 2 . 
We use P 2 /R 0 = O.71(M2/M0)°-^^ for a star in thermal equi¬ 
librium also from HNR. These M 2 -R 2 relations give a pe¬ 
riod gap (transition from the former to the latter) at 2-3 hours 
when the donor is 0.2 IM©, as found for this “standard model” 
(Model A) in HNR. 

In addition to Mign and the IMB prescriptions we also re¬ 
quire a CV population model to construct the CN orbital pe¬ 
riod distribution. Detailed Monte Carlo CV populatio n mod¬ 
els are explored by HNR and by iNelson et al.l (120^), in the 
same context as here, but without the TB04 predictions of U 
from M. HNR calculated the number of CVs per period in¬ 
terval, which we denote np so that npdP is the number with 
orbital periods between P and P + dP. We use a single char¬ 
acteristic mass in lieu of a detailed Mwd distribution. For 
the HNR population we use Mwd = 0.72M©, the average CN 
Mwd for their populations, and Z)He = 0.001. 

Using the above relationships to obtain M, np from HNR, 
and the Mign from TB04, the resulting CN distribution is 
shown by the dashed line in Figure |2] and is normalized to 
match at Porb = 3.5 hours. The most noticeable difference be¬ 
tween this and our steady-state model presented below (thick 
solid line) is how the CV population based on HNR drops off 
at long Porb due to the CVs being born throughout the 3.5- 
6 hour range. The exact fo rm of this turnover depends on 
the common envelope phase JNelson et alJl2004l) through the 
orbital period distribution of the common envelope products. 
By construction, the HNR distribution lacks systems in the 
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period gap. Every system which comes into contact with a 
companion that has a radiative core (Eorb ^2.4 hours) goes 
through a MB phase. This distribution matches quite well 
in both the 1-2 hour and 3-4 hour regions, keeping in mind 
the ^/N noise present in the observed distribution. The HNR 
population extends down to 1 hour due to the famous discrep- 
ancy in the period min imum between observation and theory 
dKolh&Raraffel 1999ft . 


3.2. Simpler Population Calculations 


In order to understand how the Porb distribution depends 
separately on the distribution of M and the CV birth distribu¬ 
tion, we also construct a steady-state population with a spec¬ 
ified birth distribution. CVs are born (i.e. come into contact 
and initiate mass transfer for the first time) at a variety of or¬ 
bital periods determined by the distribution of M^o and M 2 
after the common envelope event. For a population of np CVs 
per period interval with characteristic mass Mwd, conserva¬ 
tion implies 


dnp _ 5(n/>Porb) , 


(3) 


where bp is the CV birth rate per period interval, i.e. the num¬ 
ber of CVs per year which commence accretion for the first 
time between P and P + dP. We assume that all matter ac¬ 
creted from the companion is ejected in the CN, consistent 
with current observations (TB04), so that Mwd is constant. 
We also assume that the CV birthrate has been constant over 
the last 4 Gyrs, the time it takes a CV to reach the period 
minimum at Po^b — 80 minutes (HNR). Beyond this point the 
companion becomes essentially fixed in radius, so that the bi¬ 
nary expands and M drops quickly. These assumptions give a 
steady-state (dnp/dt = 0) CV population whose number den¬ 
sity per period interval from equation (□ is 

2 pPoTh 

np = - / bpdP, (4) 

I-* orb I J cxD 

where we assume no systems at large Porb- 
As noted above, the secondary’s mass-radius relation is 
well characterized by a power law, and so for Roche Lobe 
filling systems a = £/lnPorb/'^lnM 2 is approximately a con¬ 
stant, typically between 0.5 and 1. Thus Porb = -MaPoih/M 2 
and the CN frequency per period interval is 


Vp = 


M 


np oc 


442/Porb 

MigYi 


(5) 


The direct dependence on M cancels for this observable. By 
far the most important part of the remaining M dependence 
is in Mign. Evaluating across the period gap, M 2 is constant, 
while Porb changes by 50% due to bloating of the companion 
above the gap. On the other hand M is expected to change 
from 10~‘^Mq yr“' to ~ 5 x 1O“"M0 yr“*, causing Mign (from 
Figure [0 to increase from ~ 1O“^M0 to > IQ~‘^Mq, a fac¬ 
tor of > 10. Thus the CN orbital period distribution directly 
reflects M through its effect on Mign. With the consistent cal¬ 
culation of Mign from M provided by TB04, we can probe the 
M-Porb relation directly by comparing the steady state popu¬ 
lation to observations. Evidence for evolution across the gap 
is then contained in the degree to which a single steady state 
distribution can reproduce the observations. 

We use a binary with Mwd = 1 -OMq as representative, since 
the actual mass distributi on is very sensitive to the com¬ 
mon envelope prescription JNelson et alJl2004l) . Also, we set 


Whe = 0.001 because this is the approximate level expected 
in the accreted matter (See section |2}. In order to account 
for the possibility of systems which do not undergo magnetic 
braking, we use three populations. The first follows the IMB 
scenario (the IMB population), while the others only have Jgr 
applied (the NMB [No Magnetic Braking] populations). In 
the NMB case the secondary mass-radius relation is that of an 
equilibrium main sequence star from HNR. 

Though the detailed form of the birth distribut ion above the 
gap i s sensitive to the CE prescription (HNR; iNelson et al.l 
we employ a constant bp between 3 and 5.5 hours for 
the IMB population. This capures the contrast in classical 
nova rate across the gap without depending explicitly on a 
particular CE prescription. CVs born with Porb ^ 2.5 hours do 
not undergo magnetic braking in the IMB scenario. We call 
this population the NMBb (NMB below gap) population, and 
it has bp constant between 1.3 and 2.5 hours. A second NMBa 
(NMB all periods) population represents systems which do 
not undergo magnetic braking regardless of their birth Poib 
(see discussion of magnetic systems below). This population 
has a constant birth function between 1.3 and 5.5 hours. Note 
that, unlike the HNR population, we directly specify the min¬ 
imum Porb. These birthrates can then be used in equation 0 
and 0, and the relative birth fractions for each population are 
given by bp{P^^^ 

Emin)' 

This steady state population is shown as the thick solid line 
in Figure|2 arbitrarily normalized to match the observed dis¬ 
tribution (thin solid line) at 3.5 hours. The relative birth frac¬ 
tions for the curve shown are NMBa; 8%, IMB; 46%, NMBb; 
46%. This model reproduces the major features of the ob¬ 
served population below 4 hours, including the slope in the 
cumulative distribution between 3 and 4 hours. The depen¬ 
dence on the common envelope prescription and binary popu¬ 
lation parameters is essentially contained in the average Mwd, 
to which the shape of the curve is insensitive. The effect of 
a 30% change in Mwd is small compared to the > 20 times 
difference in M across the period gap. 

The contrast in CN rate across the period gap does have 
some dependence on VsHe- Under gravitational braking only, 
as is the case below the period gap, M is low enough that Mign 
is insensitive to Wjje- However, the CN rate above the pe¬ 
riod gap, where magnetic braking enhances M, does depend 
on Wjje. The dependence can be estimated by using the ig¬ 
nition density, pign oc from TB04, and combining 

with an ideal gas EOS, Mign cx Eign oc PignTign, and the power- 
law form of a free-free radiative envelope T oc jq 

tain Vp IX. I /Mign oc This scaling is confirmed by nu¬ 

merical calculations. Thus the expected maximum value of 
Wjje = 0.003 (see Section|2} leads to a 50% enhancement of 
the CN rate above the period gap over Vajje = 0.001. Matching 
this new normal distribution to the observations at 3.5 hours, 
the rate below the gap is 2/3 of that shown in figure|2] still con¬ 
sistent with the observations. Thus the value we have adopted, 
Wjjg = 0.001 is expected from donor evolution and provides 
the best fit to the data when birthrates above and below the 
period gap are assumed to be approximately equal. 

The observed orbital period distribution of CN is consistent 
with CVs evolving across the period gap, and magnetic brak¬ 
ing between 3 and 4 hours governed by a law like equation^] 
Ai shown by the thick (model) and thin (data) solid lines in 
figure^ the contrast between the 1-2 and the 3-4 hour ranges 
is of the appropriate size for a population in which J above 3 
hours is enhanced by magnetic braking. This contrast in CN 
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rate is entirely due to the dependence ofMig^ on M. 

3.3. Magnetic Accretors with No Magnetic Braking 

The two components included in the model population, 
IMB and NMB, are also shown separately as the dot-dashed 
and dotted lines. These demonstrate the significant enhance¬ 
ment that magnetic braking implies in the CN rate over gravi¬ 
tational radiation alone. The prediction of a NMB model with 
bp = 0 (all systems being born at long periods) falls weakly 
to shorter orbital period in the normal distribution, as for the 
IMB curve below 2 hours. With a constant bp, a NMB model 
monotonically increases to shorter orbital periods as seen for 
the dotted lines. Thus a NMB model alone is inconsistent 
with the observed orbital period distribution. The period dis¬ 
tribution is also satisfactorily reproduced without the addition 
of the NMB population, making our confirmation of the IMB 
prescriptions and evolution across the gap robust with some 
evidence for additional systems in and below the gap. Mea¬ 
surement of this additional contribution is limited by the cur¬ 
rent quality of data. 

There are two reasons a CV might not undergo the MB 
phase. For example, if MB cannot work when the companion 
is fully convective, then there is a minimum companion mass 
(O. 25 M 0 ) for having an MB phase. We have implemented this 
by placing all systems born with Po^b < 2.5 hours in a NMB 
population. This can also be seen in the the small number of 
systems between 2 and 2.5 hours in the HNR population. Ad¬ 
ditionally, r notivated by th e lack of a clear p eriod gap for mag¬ 
netic CVs (IVerhiinrill 997h . iTi et al.l (l1994h have conjectured 
that in Polars the WD magnetic field prevents MB. Recent 
Teff rneasurements of magnetic WD s (lArauio-Betancor et al.l 
l20()5t iGansicke & Townslevll200-^ show that above the pe¬ 
riod gap, magnetic CVs have a significantly lower M than 
Dwarf Novae at the same orbital period, providing direct evi¬ 
dence that magnetic CV evolution is quite different from that 
of non-magnetic systems. 

The absolute relative populations of magnetic and non¬ 
magnetic CVs is not well known because they are dom¬ 
inated by different and difficult to quantify selection ef- 
fects. Approximately 22% of known CVs are magnetic 
llArauio-Betancor et al.ll200^ . Using this steady state two 
component population model we find that, due to the slower 
evolution in the absence of MB, this requires that 8% of CVs 
born have strongly magnetic WDs, comparable to the fraction 
among field WDs (lTordanll997t [Wickramasinghe & Ferrari^ 
I 2 OOOI) . While giving tantalizing hints, the small number of ob¬ 
served systems makes it impossible to separate the prospec¬ 
tive magnetic and non-magnetic contributions to the NMB 
population. 

4. DISCUSSION AND CONCLUSION 

We have shown that a very simple CV population model 
with the assumptions of IMB and the Mign calculated by TB04 
reproduces the orbital period distribution of CN, supporting 
the idea that CVs evolve across the period gap. This agree¬ 
ment is independent of the outcome of the common enve¬ 
lope phase or properties of the primordial binary population. 
CN are one CV phenomena that is sensitive to M averaged 
over the longest timescales, the CN recurrence time, 10^-10’ 
years. The agreement with IMB provides additional evidence 
that the IMB prescriptions do indeed provide accurate long 
term average M’s for CVs. iTownslev & BildstenI (l2003h sim¬ 
ilarly found that the WD Teff in dwarf nova systems, which 


reflects M due to compression of the envelope, matches the 
predictions of IMB. The existence of CN with Porb in the 
period gap implies that not all systems undergo IMB evolu¬ 
tion, possibly due to inhibiti on of the seco ndary’s wind by 
a strong WD magnetic field (ITi etalJI1994l) . This inhibition 
will lead to slower evolution of magnetic CVs, and thus in¬ 
crease their relative number in the CV population. We find 
that in order to obtain the 22% fract ion of magnetic CVs ob¬ 
served llArauio-Betancor et al.ll200^ . a relative birth fraction 
of magnetic CVs o f 8% is required, comparable to the fraction 
among field WDs (lTordanll997HWickramasinghe & Ferrariol 

IMO). 

With a satisfactory M distribution in hand, the total number 
of classical novae observed provides a measurement of the 
CV population. This conversion does not depend directly on 
CE parameters, because a relatively small number of both ob¬ 
served novae and CVs in general have Porb > 4 hours, where 
the birth distribution is important. The conversion does, how¬ 
ever, depend on the typical WD mass, due to its importance 
in setting Mig„. Note also that while most CV WDs have 
Mwd ^ 0.7Mq, including many helium WDs (HNR), the av¬ 
erage CN mass is sensitive to the upper end of the mass distri¬ 
bution dNelson et alJl2004l) . Knowledge of the M distribution 
allows proper adjustment for the fact that while most CVs are 
below the period gap, most CN are in CVs above the period 
gap. This is a fairly short-lived period of the CV’s lifetime, 
lasting only ~ 100 Myr (HNR), which is preceded by the post 
common envelope inspiral and followed by accretion below 
the period gap during which the CN rate is much lower. We 
find that one CN per year implies 3 x 10^ (9 x 10^) CVs above 
the period minimum, t <3 Gyr (4 Gyr) and a CV birthrate of 
10“^ yr“' (2 X 10“"^ yr^) if the average CN mass is I.OMq 
(0.6Mq). 

Since pre-CVs finish the common envelope phase with 
a wide range of orbital periods, the total CV and pre-CV 
population, and thus the CN rate, should relate to the to¬ 
tal amount of low-mass stars, which is well represented by 
the K-band light. The apparent dependence of nova rate 
on stellar population is still uncertain, but points to these 
old er populations. Surveys o f nova rates in external galax¬ 
ies dWilliams & Shaft^l2r)04li continue to show that the K- 
band specific nova rate is about 2 ± 1 x 10“'** Lqk 
does not depend strongly on galaxy type in nearby galaxies, 
with the exception of the EMC and SMC. M31 shows evi¬ 
dence that the cla ssical nova rates migh t be dominated by the 
bulge population dShafter & Irbvll2()()ll) . Our work allows us 
to convert the K-band specific Nova rate into the CV pop¬ 
ulation and birthrate, giving 60-180 CVs above the period 
minimum for every 10°Lq k in an old stellar population and 
2(4) X 10-4 CVs born per year per lO'^Lo,/^. Th e 24-133 CVs 
estim ated from Chandra observations of 47 Tuc dHeinke et al.l 
l200.5h is consistent with t he 60-180 predicted by our work for 
the cluster mass (10®Mq iPrvor & Mevlanin~993l) . and agrees 
with the CN rate implied by the one or two CNe observe d 
in the last ~ 140 yr in the Galactic GCs dShara et al .1120041) . 
Thus the GC CV population appears comparable to that in 
the field. The Type la supernova rate in elliptical galaxies i s 
3 . 54 J [ X 10-4 per year per IO'^Lq^a- dMannucci et alJl2004l) . 
remarkably close to the CV birthrate. However, all indications 
are that the WD masses in CVs are not nearly large enough 
(TB03) to be Type la progenitors dBranch et al 

The loc al galactic K-band lumin osity density is 0.1-0.2 
Lq^k pc-^ dolling & Merrifieldl2001h . giving a local CV space 
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density of approximately 9(27) x 10“^ pc“^ again for an av¬ 
erage mass of 1.0 (0.6 ) Mi^,. This ap ees well with the 

2 X 10“^ pc“^ found by iPolitanol (119961) from a theoretical 
calculation of the CV birthrate. This space density is only 
slightly higher than the 6 x 10“^ pc“^ inferred from the PG 
survey lRingwal dll996l). but mo re than 10 times the 8 x 10“^ 
pc“^ obtained bv l""e idUH). Both these survey numbers 
assume a scale height of CVs of 150 pc, whereas the above K- 
band luminosity density assumes a scale height of low mass 
stars of 300 pc. Although more modern color selected surveys 
exist, the lack of a broadly applicable, accurate CV distance 
measure hinders improved measurements of the space density. 

In a similar manner, the total rate of matter injected into the 
ISM can be evaluated in terms of the CN outburst rate. Since 
we are assuming that all the matter transferred onto the WD is 
ejected in each CN outburst, this is just the average M for the 
population times the total number of systems. We find that 

3 X 1 O“^M 0 yr“' (9 x IO'^Mq yr“') is ejected into the ISM 
for each CN/yr if the average CN mass is 1.0 Mq (O. 6 M 0 ). 
The average initial mass of the donor contributes additional 
uncertainty, but generally less than that due to the unknown 
average WD mass. Though somewhat lower than the value o f 
2 x 1 O“"^M 0 per CN outburst found bv lTriiran & Tiviolill986l) . 
these values are sim ilar to the ~ 1 O“'^M 0 per outburst used by 
iGehrz et"^ (Il998h to estimate the impact of nucleosynthesis 
in CN outbursts on both the ISM and primordial solar system 
material. An interesting effect of the contrast in M across the 
gap is that while approximately 85% of CN outbursts occur 
above the period gap and only 2 % of pre-period minimum 
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systems are above the gap, approximately equal processed 
matter is contributed by objects above and below the period 
gap before the period minimum. This means that the type of 
CN which are observed most often are only half the story for 
ISM injection. 

With proper calibration, classical novae provide a use- 
ful standard candle for ex tragalactic distance measurements 
ddella Valle & Livioll 199.51) . The principle impediment to this 
usage is that the single measured parameter available for cali¬ 
bration, the rate of decline, depends on Mwd, M and Tc- TB04 
calculated Tc consistently from M, eliminating it as a free 
parameter. Here we have demonstrated that IMB provides 
a satisfactory prediction of the M distribution in the Galac¬ 
tic CV population, finding that CN outbursts occur predomi¬ 
nantly when M ^ yr“', with a significant tail at lower 

M. This leaves only Mwd to be calibrated from the rate of 
decline. Thus, with the support of theory, a statistically com¬ 
plete measurement of nova magnitudes and decline rates may 
be sufficient to measure the distance to an external galaxy. 
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